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ABSTRACT 



Coherent light beams propagating through the atmosphere undergo 
considerable phase perturbations due to fluctuating temperature structures in 
the atmosphere. Understanding and measuring these structures on a real-time 
altitude dependent basis is inherent to successful deployment of ground based 
lasers and particle beams. 

One method used to detect these temperature structures is an acoustic 
profiler, or echosounder. Of immediate interest is the ability of high frequency (5 
kHz) planar array echosounders to rapidly detect low level turbulence (below 
200 meters) and quantify the results. 

This thesis involves design improvements of previously developed 
echosounder arrays and associated software. Particulary, this thesis 
demonstrates that tighter packing of elements in a planar acoustic array 
produces better side lobe reduction than less densely packed arrays. This 
results in higher energy density in the main lobe and increased performance. 

Also included in this thesis is a computer method which allows relatively 
accurate beam pattern prediction from any given planar array. 
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I. INTRODUCTION 



Optical degradation of coherent light in the atmosphere is strongly a 
function of small scale temperature gradients which change the local index of 
refraction in the atmosphere [Ref. 1]. These localized structures in the 
atmosphere fluctuate rapidly, changing in size, density, location, velocity and 
interaction with neighboring structures. There is also considerable variability 
with diurnal temperature changes, annual weather variations and regional 
climatological differences. Understanding the mechanisms of this turbulence 
and measuring it on a real-time local level are intrinsic to high power earth 
based lasers. 

These small scale turbulence structures in the atmosphere play an even 
greater role with respect to acoustic absorption and reflection [Ref. 2]. It is well 
known that acoustic pulses, on the order of 1 to 6 kHz frequency, are ideal in 
defining turbulence structures in the first few kilometers of the atmosphere, with 
lower frequencies having greater ranges. Until recently, the most popular type 
of acoustic sounder was a monopole speaker/receiver mounted in a parabolic 
dish [Refs. 3 and 4], Several problems with this method include difficulty of 
transportation, high amplitude side lobes at large angles and poor altitude 
resolution due to the low frequencies involved. These deficiencies can be 
overcome by means of planar arrays [Refs. 5 and 6], and this thesis was 
designed to amplify this point. 
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Working under Walters, Weingartner and Wroblewski [Refs. 7 and 8] again 
demonstrated the feasibility of this technique by building and testing a 25 
element 5 kHz planar array. The result was an extremely effective low altitude 
echosounder which was less expensive, more portable and had better side 
lobe reduction than previous monople sounders. 

During the course of this previous work several questions arose 
concerning the planar array. If the elements were more densely packed, then 
could the size of the side lobes be significantly reduced, thus increasing the 
accuracy of the beam? Could the array be modeled on a computer to produce 
the best design? Could the enclosure be improved? 

This thesis is an attempt to answer these questions, and involves the 
development of a dense packed 5 kHz hexagonal array and associated sound 
deadened enclosure. Additionally, various computer programs using fast 
Fourier transform algorithms were studied as a first approximation in 
diagnosing acoustic side lobe suppression. 

The results indicate that densely packed planar acoustic echosounders 
produce better side lobe suppression than more open arrays and are accurate 
tools in the identification of low level turbulence. 
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II. BACKGROUND 



The atmosphere induces random phase perturbations on an 
optical beam. Diffraction of these small scale phase variations produces a 
larger beam, restricting the ability to focus energy on a target. Optical phase 
perturbations in the atmosphere stem from temperature fluctuations carried by 
the turbulent velocity field. These temperature fluctuations vary in size, 
cascading from tens of meters in diameter down to centimeters and smaller. 
These structures appear to be most variable and inhomogeneous at the 
interface between stratified layers, particularly near the earth's surface. 

Within the first few hundred meters of the surface, these temperature 
structures are dominated by the surface heat flux, resulting from differing rates 
of heat transfer between the earth and the atmosphere. During the day the 
surface heats much more rapidly and non-uniformly than the atmosphere. An 
insulating blanket of air near the surface heats and rises off the ground, more or 
less in sheets. Warmer sections rise more rapidly, quickly disrupting the 
uniformity of these sheets. Cooler air moves in to replace the rising air, further 
disrupting the sheets and causing vertical thermal plumes to form. These 
thermal plumes continue upwards in the atmosphere, slowly mixing and 
increasing the temperature of the lower atmosphere. This cycle is continuous, 
with the actual transfer of energy occurring by small scale molecular diffusion. 
[Ref. 9]. 
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At night the cycle reverses, with the surface cooling off much more rapidly 
than the air via Planck radiation. In this case, as the air gives up energy to the 
cooler surface, it produces a turbulent boundary layer (via velocity interaction 
with a stationary surface) which remains relatively close to the ground. This 
cooling process is generally less turbulent than the heating process, due to the 
reduced heat flux at the boundary, although intense turbulence continues to 
exist at the interface of the horizontally stratified layers. 

There is also a neutral event, during the transition between night and day, 
where the surface temperature equals the ambient air temperature. During this 
time the net heat flux between the air and the surface changes sign and passes 
through a minimum. The atmosphere generally becomes optically stable during 
these brief periods. 

The structural variability is clearly evident in optical instruments designed 
to measure the turbulence of the atmosphere [Refs. 10 and 11], but altitude 
correlation with optical instruments alone is difficult. Acoustic echosounders 
alleviate this difficulty. 

Acoustic echosounders are devices which generate and propagate 
energy into the atmosphere and measure the backscatter return from turbulent 
structures as a function of altitude. The power returned from these turbulent 
structures, Pr, is proportional to the power transmitted, Pj. The equation 
relating these terms has been summarized by Neff [Ref. 12], based upon work 
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by Tatarski [Ref 1], and is repeated in equation (1) below. 

(1) P R = P T E r E t (ct/2) {exp(-2IBR)} (AG/R 2 ) p[R,/] where, 

- P R is the electrical power returned from range R, 

- P T is the electrical power transmitted at frequency /, 

- Er is the efficiency of conversion from acoustic to electric power, 

- E T is the efficiency of conversion from electric to acoustic power, 

- c is the local speed of sound (m/s), 

- 1 is the acoustic pulse length (s), 

- A is the antenna area (m 2 ) , 

- R is the range to the scattering volume (m) 

- exp (-2I3R) is the round trip power loss due to attenuation where 13 is the 
average attenuation to and from the scattering volume, 

- G is the effective aperture factor for the antenna, 

- p [r,/] is the scattering cross-section per unit volume; that is, the 
fraction of incident power backscattered per unit distance into unit solid 
angle at frequency /. 



The backscattering cross section, p, is the unknown quantity of interest in 
Eqn. (1). Tatarski [Ref. 1] uses p in equation (2) (below) to find Cj2, the 
temperature structure parameter. 

(2) p[R,/] = 0.0039 k 1/3 (Cj 2 / T 2 ) where, 

- k is the acoustic wave number = 27r / L, 

- L is the acoustic wavelength, 

- T is the absolute temperature and 

- C T 2 is the temperature structure parameter. 



The temperature structure parameter, C T 2, is the mean squared 
temperature difference between separate points in space, 

(3) Ct2 = {<T 2 (R2)-T 1 (R 1 )>2}/(R 2 -R 1 )^ , 

and provides the scattering centers which produce backscattered energy from 
the turbulence. <AT> 2 empirically increases by AR 2/3 . Therefore, it is normalized 
by AR 273 to produce a relatively constant temperature structure parameter. 

Since p is proportional to the backscattered return power, P R , and to 
the temperature structure parameter, C T 2, then P R is proportional to Cj2. 
Solving equation (1) for p and equating equations (1) and (2), the following 
expression for Cj 2 is given for an acoustic echosounder of wavelength L; 

(4) C t 2 = 1/0.0039 (1/ E r E t ) T 2 k‘ 1/3 (2/ct)(1/AG) P R /P T {R 2 exp(2f3R)}. 

The electric to acoustic efficiency factors for the transducers (speakers) 
were measured in Reference 7 and are inherent to the design of the speakers. 
E R = Ej= 0.5, an excellent value for speaker efficiency, which accounts for 
their use. The aperture area, A, for the 19 element hexagonal array is 0.0866 
m 2 . 

The effective aperture factor for the array, G, accounts for the non-uniform 
antenna directivity. G is a difficult quantity to calculate or measure. It represents 
a normalized average of the acoustic energy within the transmit/receive 
antenna beam lobe and is clearly summarized by Probert-Jones [Ref. 13]. 



Based upon Reference 13 and studies by Hall and Wescott [Ref. 14], a value of 
G = 0.4 is currently used in computations. This quantity is presently undergoing 
extensive research by Dr. Walters. 

Another difficult factor in Equation (4) involves the average atmospheric 
attenuation loss factor, 13, which is strpngly dependent upon frequency and the 
absolute water vapor in the atmosphere [Refs. 2 and 15]. This dependence is 
illustrated in Figure 1, which is taken directly from Reference 15. The water 
vapor pressure is measured via a digital humidity/temperature indicator, model 
5165-A from Weathermeasure [Appendix A], 



WRTER-VRPOR PRESSURE mb 

Fig. 1. Attenuation versus Water-vapor Pressure (mb.) 
for Frequencies around L.6 KHz. 




Combining these known values into a constant Q, equation (4) becomes; 

(5) CyS = (Q) T2k*i/3 {P R / (ct)( P T )} R 2 exp(2(3R) where, 

Q = {2Pr}/{0.0039 P T Er Ej A G}. 

Wroblewski [Ref. 8] utilizes this in the software for the echosounder to 
calculate a 15 minute time averaged value for the temperature structure 
parameter as a function of altitude. Uncertainties in the antenna beam shape, 
side lobes, rejection, electronic gains, transducer efficiencies, electronic 
antenna gains and the atmospheric attenuation produce uncertainties in the 
absolute calibration of an echosounder. 

Independent verification of the C T 2 values measured and calculated by 
the echosounder are necessary to confirm the accuracy of this technique. 
Without this confirmation it is still possible to guage relative variations in the low 
level turbulence, simply with the echosounder profiles. These profiles are 
presented and discussed in the following chapters. 



III. EQUIPMENT 



Much of the equipment used in this project was borrowed from 
the 25 element square array designed by Weingartner and 
Wroblewski. The equipment basically involves a Hewlett Packard 
(HP) 217 computer for data generation and analysis, an HP 331 4A 
function generator (controlled by the software) which produces the 
initial pulse, an amplifier, a switch/pre-amp which passes the 
initial pulse then amplifies and transfers the return signal to the 
computer, the array and shroud which act as the physical 
transmitter/receiver and a filter used to eliminate stray noise. 
These units will be broken down and individually discussed in this 
section. The software developed for the echosounder by Wroblewski 
was not significantly altered and is not included in this document. 
[Refs. 7 and 8] 

A. ELEMENTS 

The elements of the array consist of 19 piezoelectric 
Motorola brand speakers, model # KSN 1005A, with attached horns 
of model # KSN 1032A. Ceramic piezoelectric drivers were chosen 



because of their small size, excellent (see Figure 2 ) frequency 
response at 5 kHz, and a shorter 'ring' time with respect to 
electrodynamic speakers. Any speaker acts as a damped 
harmonic oscillator with a 'ring' time (recovery time) inversely 
proportional to the damping constant. The damping is much stronger 
in Piezoelectric speakers than in conventional electrodynamic 
elements. 
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Figure 2. Speaker Frequency Response 



An initial set of 40 speakers was tested to find 19 which 
were similarly matched in frequency response. This testing 
was done by comparing the output of a chosen reference 
speaker with the acoustic-to-electrical conversion of each 
test speaker. By clamping the reference speaker face-to-face 
with the test speaker, a known voltage to the reference should 
produce a range of voltages from the test specimens. These voltages 
were compared using Lissajous plots on a Nicolet Oscilloscope. It 
was more important to obtain a set of identical speakers rather than 
the best speakers, in order to achieve a more uniform final array. 

The horns on each driver, originally circular with a square 
flange face were trimmed to a circular diameter of 7.62cm such that 
they fitted close packed. Ideally a diameter of d= L/2 
(L=wavelength) would provide optimal phase spacing of elements, 
but this required either smaller speakers ( d = 3.4 cm ) which 
were unavailable, or a lower frequency (=2230 Hz), which was 
below the useable range of available speakers. 

Another source of speaker was identified as a Philips brand 
piezoelectric 2002 PT. [Ref. 16] with better frequency response at 



3kHz than the Motorola speakers In use. Unfortunately these Philips 
speakers are no longer manufactured. 

B. ARRAY 

The speaker elements were configured in a close-packed 
planar array on an aluminum-covered balsa wood board to absorb 
vibrations (Fig. 3). The elements were wired in parallel, reducing 
the array impedance to 300 Cl / 19 = 15.8 Cl. The mounting box was 
locally constructed of sheet aluminum and enclosed a layer of 
semi-porous foam loaded with a thin lead sheet to reduce stray 
noise signals. Input/Output was via a coaxial connector. Bracketts 
were attached to allow 360° rotational testing from 2 axes in the 
anechoic chamber. 




Figure 3. 



Array Dimensions 



C. ENCLOSURE 



The enclosure was designed to be easily reproducible, 
inexpensive, portable, and easily assembled [the earlier enclosure 
took some time to set up and break down]. Optimally, a long 
cylindrical enclosure was desired to reduce side lobe scatter and 
return without obstructung the main lobe. The enclosure chosen to 
optimize the above requirements was a Rubbermaid 55 gallon plastic 
trash container, with a height of 84 cm, having a slight upward 
taper; the barrel has a lower inside diameter of 53 cm, opening to a 
diameter of 60 cm. This provided a finished open half angle of about 
20°, 8° more than optimally desired according to tests in the 
anechoic chamber with the unencumbered array [see RESULTS'! . A 
second possibility was a highway barrier barrel [Appendix A], but 
these were much more expensive and involved similar dimension 
problems. 

The inside surface of the enclosure was covered with lead lined 
semi porous foam, with 7.62 cm convoluted (egg crate shaped) foam 
on one side, and 2 cm flat foam on the backing side. The entire foam 



liner was attached with standard contact cement. The lead lining 
amounted to a density of about 1 Ib/sqft which is sufficient to 
reduce one way noise penetration by 40dB. The convoluted foam 
helped break up internal and off axis scatter (spectral reflection) 
and aided in side lobe reduction. The foam used turned out to be the 

last of a salvage batch of foam, and a second source of foam was 

desired. An ideal substitute for lead lined foam (which was 
difficult to obtain) was found to be a barium impregnated vinyl foam 
(a thin rubbery vinyl layer loaded with barium and sandwiched 
between semi-porous foam) with a similar density of 1 Ib/sqft and 
virtually identical sound absorption capability (40dB). Two 
sources of this foam were identified [Appendix A], 

D. ELECTRONIC HARDWARE 

The data processing system consisted of an HP Model #217 
computer with a 20 megabyte hard drive and 2 megabyte memory, an 
INFOTEK BC 203 Basic compiler, analog to digital (A to D) converter, 

INFOTEK FP 210 floating point accelerator (to enhance the speed of 

execution), HP 9133 floppy disk drive, monitor and HP ink-jet 



printer (see Figure 4). The HP 217 was the workhorse of the system, 
triggering the function generator, receiving and processing the 
return pulse, and outputting the results in graphical mode. Without 
this unit, timely acquisition and reduction of data would be 
impossible. 

The HP 3314A function generator, adjusted and triggered 
automatically by the HP 217 computer, was used to generate a 
specific transmitted signal. This signal was generally 100 cycles of 
a sinusodial waveform of 1.5 volts (rms) amplitude. 

A QSC model 1700 audio amplifier was used to increase the 
voltage to the array from 1.5 volts (out of HP 3314A) to 30 volts 
(rms). This model was chosen because of its ability to drive a 
highly reactive load (the array). A lighter, more compact amplifier 
could easily replace this item of equipment, and will be considered 
during future upgrades. 

A pre-amplifier (pre-amp) (Fig. 5) was used to isolate the 
transmitted pulse from the received pulse and to amplify the return 
signal by a gain of approximately 11,000. The pre-amplifer was 
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TO AND FROM ARRAY 



RETURN SIGNAL (TO COMPUTER) 




FOR FILTER, USE 
INPUT GAIN 1 AND 
OUTPUT GAIN 10. 



Fig. 4 Echosounder System Set-up 
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TRIGGER SIGNAL 




Fig. 5. Schematic of Pre-Amplifier 
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designed and built by Walters and utilizes two low noise solid state 
operational amplifiers, an LT 1007 and LT 1037. All four isolation 
rectifiers were of type IN 4000. 

The important properties to consider for an effective pre-amp 
were complete isolation of the circuit during transmission (to 
prevent overloading components) and immediate signal retrieval 
following the 'ring' time of the array. The more rapid the recovery 

time of the array (a function of driver capacitance (C a ) and the loop 

resistance (R i )) , the quicker the pre-amp can respond to returned 
signals and the smaller the surface blind zone. For this array, the 
capacitance of each driver was Cq = 0.1 j i F. Thus the time constant 
for the array was 

RC = (C[))(#drivers)(Ri ) = time constant 
= (0.1*10- 6 F)(1 9) (270 Q) = 0.513 msec. 

For a given array, resistor R-| provides the time constant of the 
circuit. 

After 5 time constants there is virtually no 'ringing' left on the 
array, and this sets the minimum height of the blind zone. Thus: 
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5 x T = 2.57 ms and H m j n = (340 m/s)(5.32*1 0'^s) = 0.87 meters. 

Therefore the blind zone stands at approximately 1 meter above the 

array, which is not noticeable on the echosounder printout. 

A Rockland Wavetek model #753A filter {115 dB/octave} was 

used to isolate noise and further amplify the return signal prior to 

computer analysis. The input side of the filter was set up as a low 

pass filter at 5 kHz and gain setting of O dB (xl); the output side 

was set up as a high pass filter at 5 kHz with a gain setting of 20 dB 

(xIO). This produced a band pass filter centered at 5 kHz with a 

band width of about 100 Hz at the *3 dB level. The over all gain of 

the filter in this mode was measured as 

output 685 m volt 

= = 9.67 = 10 

input 70.88 m volt 

The combination of the pre-amp and the filter produced an overall 
return signal gain of (pre-amp gain)(filter gain) = (1 1 094)(9.67) 

= 107,000. 

The gain of the system is thus: 107,000 * 270 Q/( 270 + 15.8 a ) 

= 101 , 000 . 
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IV. RESULTS 



A. LOBE PATTERNS 

The array was tested in the anechoic chamber utilizing a 
program designed by Butler [Ref. 17]. This program utilizes the 
set-up shown in Figure 6 and allows the measurement of beam 
patterns for various arrays. A band pass filter was used to increase 
gain and reduce 60 Hz noise from the equipment. Present on all the 
lobe patterns is a spurious lobe at * 135°. This was determined to 
be a reflection from the inside door handle of the anechoic chamber 
and should be ignored. All lobe patterns have been normalized to the 
maximum voltage detected, allowing easy comparison between 
figures. 

The hexagonal array was rotated around two axes to obtain the 
maximum sidelobe field of the array (Figs. 7 and 8). The pattern for 
the square array is shown for only one axis of orientation (Fig. 9). 
Comparison between the two arrays (Figs. 7 and 9) shows a definite 
decrease in side lobe number and intensity for the hexagonal array 
versus the square array. 
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ROTATOR 




ROTATOR CONTROLLER: BATT IS "OFF"; CCW/"CW", ON/W" 
FULL SCALE/ZERO / "OPERATE.” 

723: SET TO READ •OHMS' < 2WIRE fi>; ADJUST THE ARRAY SO 
THIS STARTS IN THE "OYLD" (OYERLOAD) POSITION 
724: SET TO READ "YOLTSAC" (~Y). 



Figure 6. Anechoic Chamber Set-up 
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19 ELEMENT HEXAGONAL ARRAY - APEX MOUNTED 




MIKE 



RUN 97 MIKE- 6.2938 V INPUT-2.5 V FREQ-5050 HZ FILTERED 



FIG. 7. HEXAGONAL ARRAY BEAM PATTERN - APEX MOUNTED 
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19 ELEMENT HEXAGONAL ARRAY - SIDE MOUNTED 




MIKE 



RUN 90 MIKE- 6.4951 V INPUT-2.5 V FREQ-4950 HZ 355 PTS FILTERED 



FIG. 8. HEXAGONAL ARRAY BEAM PATTERN - SIDE MOUNTED 
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5X5 ACOUSTIC ARRAY BEAM PATTERN 




R=3. 60 N MIKE VT= 8.2465 INPUT=*5. 0 V FREQ=5. 0 KHZ 



FIG. 9. SQUARE ARRAY BEAM PATTERN - SIDE MOUNTED 
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While the array was being tested in the anechoic chamber a 
demonstration of its detecting ability was desired. Using a heat 
gun, the array was set up in the chamber exactly as it would be in 
the field, except a Nicolet model #3091 digital Oscilloscope was 
used to detect the return signal. 

The heat gun was tested at three different settings; first was 
'heat and air', then 'air only' and last 'heat only'. With the heat and 
air both on, a definite fluctuating signal could be detected on the 
oscilloscope screen. With 'air only' there was virtually no 
detectable signal and with 'heat only’ the signal was almost back to 
its level of 'heat and air’. No measurements were taken, but the 
demonstration was repeated several times. This directly illustrated 
that the turbulence structures in the atmosphere are mainly a result 
of temperature fluctuations and not due to turbulent velocity 
fluctuations. 

The next step involved testing the side lobe suppression of the 
enshrouded array. The mechanics of this operation were 
considerably more involved than simply suspending the array from a 
rotating motor. Being too heavy to suspend, the shroud was taped to 
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a rotating turntable and centered beneath the rotator motor. The 
motor was connected to the shroud via aluminum rods and clamps. 
When the rotating motor turned, so did the shroud, allowing the 
enshrouded beam pattern to be measured and recorded with the same 
software as used previously. 

The initial beam pattern, (Fig. 10) showed that the shroud, while 
significantly reducing side lobes, was not long enough to completely 
eliminate them. Some of the side lobes were getting through. 

A first-cut step in further reducing these side lobes involved 
placing a 15 cm wide sleeve of barium impregnated smooth foam 
around the upper inside circumference of the shroud. The result (Fig. 
11) was not promising, and demonstrated that further refinement 
was required. One thing Figure 11 demonstrated, however, was that 
specular reflection from a smooth surface increases side lobe 
activity and thus convoluted foam shrouds (vice smooth foam 
shrouds) are required. 

The next step involved placing a 15 cm wide annular ring coated 
with sound deadening foam over the exit of the shroud (similar to a 
trash can lid with the center removed). This left an open circle of 
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ENCLOSED - HEXAGONAL ARRAY - FILTERED 




MIC 



RUN 115 MIC =» . 157885 V INPUT-2. 5 V FREQ=5000 HZ 



FIG 10. ENSHROUDED ARRAY 
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ENCLOSED - HEXAGONAL ARRAY - FILTERED 




RUN MODI MIC ■ . 114-418 V INPUT-2. 5 V FREQ-5000 HZ 



FIG. 11. ENSHROUDED ARRAY WITH RING LINER 
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48.25 cm diameter for the acoustic waves to propagate. This 
adjustment slightly decreased the side lobes (= 3 dB), but not 
significantly (Fig. 12). Further narrowing of the neck via more 
annular rings only succeded in reducing the main lobe, with 
little subsequent change in side lobe activity. 



B. FIELD MEASUREMENTS 

Field measurements were conducted on several occassions 
over a six month period at various locations around California. They 
generally coincided with optical measurements being conducted by 
Walters [Refs. 9 and 10] and were used to supplement the optical 
data. Of greatest interest were mountaintops above the marine 
boundary layer (=4000 ft.), where the turbulence should be a 
minimum. 

Two mountaintops in particular were visited, one being Lick 
Observatory on Mt. Hamilton and the other being Mt. Wilson 
Observatory at Mt. Wilson. These were chosen because of the 
availability of telescope domes, which were required to reduce wind 
vibration on the optical components and protect the computers from 
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FIG. 12. ENSHROUDED ARRAY WITH ANNULAR MUZZLE 
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dust and debris. The echosounder itself is impervious to wind, and 
can be operated wherever electricity is available. 

Lick Observatory was visited on 1 and 2 August 1987 and some 
of the results are presented in the following figures. All times are 
given as universal time. Subtract seven hours to obtain local 
(daylight savings) time at the site. Both days were similar with 
respect to turbulence. Therefore, only profiles from 2 August will 
be presented to illustrate the ability of the echosounder to detect 
and present the rapidly varying nature of the atmosphere accurately. 
The Cy2 profiles are currently unverified and should be ignored. 

Figure 13 was taken early in the evening and demonstrates the 
presence of convective pluming. This pluming vanishes gradually, 
eventually yielding to the nighttime structure which starts in Figure 
14. By this time a definite layer has set in and surface pluming is 
no longer dominant. Later in the evening a strong turbulence layer 
descendes (Fig. 15) which quickly settles [Fig. 16] and is quite dense. 
This layer remains close to the ground the rest of the evening, with 
light turbulence detectable above it. 
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Fig. 13. Echosounder Profile; 1830- 1845 Local Time 
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Fig. 14. Echosounder Profile; 2000 -2015 Local Time 
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Fig. 15. Echosounder Profile; 2145-2200 Local Time 
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Fig. 16. Echosounder Profile; 2215-2230 Local Time 
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Early in the morning (Fig. 17) another strong layer is noted. As 
the sun begins to rise (Fig. 18), commencement of the daily 
convective pluming becomes clearly evident. These figures 
illustrate the dramatic ability of the high frequency echosounder to 
profile low level turbulence in the atmosphere. Interestingly, 
whenever a turbulence structure was detected on the profiler, there 
was also noticeable degradation to both the optical instruments in 
use. 

The next occasion to use the echosounder occurred on 6 and 7 
November 1987 with a field trip to Mt. Wilson Observatory. 
Unfortunatley, no echosounder data could be obtained due to severe 
radio frequency [RF] interference from approximately 10 radio/TV 
transmitters co-located at the observatory. The trip was not a loss, 
however, as it was used to correct a subtle flaw in the design of the 
pre-amp. 

The problem was traced to the 'floating point ground* design 
of the pre-amp. The coaxial cables leading to the pre-amp made 
ideal antennas at the 200 megahertz range, transferring enough 
radio frequency voltage into the pre-amp casing to influence the 
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Fig. 17. Echosounder Profile; 0430-0445 Local Time 
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• Fig. 18. Echosounder Profile; 0530-0545 Local Time 
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circuit. The rectified voltage effectively overloaded the first stage 
amplifier, overriding the weak return signal from the array and 
destroying the data. This effect went undetected in previous use 
of the pre-amp. The problem was corrected by improving the 
grounding structure of the device. 

Note that all these profiles were taken with the hexagonal 
array mounted in the square shroud, as the hexagonal shroud was not 
yet completed. To ensure that the quality of data was not 
compromised, the square array was also run at the site and the 
profiles were identical. 
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V. ANTENNA BEAM PATTERN MODEL 



This thesis also modeled the planar acoustic array as an 
optical aperture, obtaining the Fourier transform of the far field 
'diffraction' pattern via computer calculations. The technique rests 
on the fact that 

. . the Fraunhofer diffraction pattern is. the Fourier transform 

of the field distribution across the aperture . . [Ref. 18, p. 412] 
The quote applies to electromagnetic radiation, but Dr. waiters 
correctly construed the electromagnetic amplitude as an acoustic 
pressure, obtaining the following equation; 

(6) P(k x ,k y ) = 1/iR JJ aperture P(x,y) exp {i (k x x + k y y)} dxdy ; 

- P(k x ,k y ) is the pressure distribution in the image plane, 

- k x and k y are the spatial frequencies along the x and y axes of 

the aperture respectively; kj= k*J/R, where J= x or y; k=2? i /\ , 

- P(x,y) is the pressure amplitude at the aperture, 

- exp {i (k x x k y y)} is the spatial variation of the aperture and 

- R is the distance to the image plane; i accounts for temporal 

phase distributions across the aperture. 
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The model can be thought of as a plane coherent acoustic wave 
passing through the aperture, which produces a diffraction pattern 
(the Fourier transform) in the far field (the image plane) by means 
of the Huygens-Fresnel principle. From this pattern it was possible 
to reasonably predict the size of the main lobe and the approximate 
angles and intensities of the side lobes. This modeling was done for 
both the square array and the more densely packed hexagonal array. 
The results are presented in this chapter. 

A. HARDWARE 

The computer used in this phase was a Hewlett Packard (HP) 
217 with an INFOTEK Basic compiler, identical to the computer used 
in the echosounder arrangement. The memory requirements for this 
program, however, were tremendous, requiring over 4 megabytes of 
memory when using the 512 x 512 grids. This lead to insertion of 2 
additional memory boards into the computer, (one from INFOTEK and 
one from EVENTIDE), bringing the hard memory to over 5 megabytes. 
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B. SOFTWARE 



The program [Appendix B] utilizes a Fast Fourier Transform 
(FFT) routine supplied by INFOTEK. The idea is to enter the aperture 
(array) as a small series of points on a huge grid. Typically the grid 
size would be 256 x 256 points square, and occassionally 512 x 512. 
The array takes up only a small center portion of this grid, 

approximately 40 x 40 grid points. The larger the grid field in 
relation to the input aperture, the better the resolution, provided the 
input aperture is not so small as to obscure the finer details. 

The aperture area was entered digitally as a series of 'ones’, 
signifying a uniform coherent plane wave at the aperture, with the 
remainder of the field entered as zeroes. The FFT routine then 

computed the far field Fourier transform of the array, and plotted 

the output intensity in relation to a gray scale. The gray scale was 
adjusted to encompass 3 decades of intensity. 

C. RESULTS AND ANALYSIS 

A sample of the input apertures, both square and 

hexagonal, are shown in Figures 19 and 20. Generally, the radius 
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FIG. 19. HEXAGONAL ARRA/ INPUT 




FIG. 20. SQUARE ARRAY INPUT 
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of each element is only 4 grid units, with a background field of 256 
x 256 grid units. These figures have been enlarged to show detail. 
The output Fast Fourier Transform for each is shown in Figures 21 
and 22 respectively. The center portion of the FFT for each array is 
the lowest spatial frequency for the array, and corresponds to 
the main lobe. Rings further out correspond to the side lobes at 

larger and larger angles. Rings far from the main lobe are not 
physical except at very high frequencies. Interestingly, as the 
frequency is increased, more and more side lobes are forced into the 
forward hemisphere. 

The program calculates the beam structure of a single piston 
driver on an infinite plane. The actual source is no! a piston driver 
(but the approximation is close). Huygens wavelettes propagating 
backwards from the end of the horns reflect from the support plane, 
increasing the energy of the actual lobe field. Also, the physical 
side lobes can extend no further than 90° from the axis of the array 
(due to surface absorption and reflection) but this limitation is 
ignored by the program. In addition, an obliquity factor was not 
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FIG. 21. HEXAGONAL ARRAY FFT OUTPUT (ENLARGED) 




FIG. 22. SQUARE ARRAY FFT OUTPUT (ENLARGED) 
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included. These limitations affect the graphical output, but do not 
significantly affect the following calculations. 

The angles of the output beam pattern rings depend on the 

frequency of interest (5 kHz), the dimensions of the actual array and 
the input/output dimension of the aperture. To calculate 

information from these printouts it is necessary to use optical 

aperture equations from Hecht and Zajac [Ref. 18]. The half angle 
formula for a circular aperture is: 

(7) sin 0 = q/R = 1 .22 L/2a where, 

- q is the distance from axis to side lobe of interest, 

- a is the radius of aperture, 

- L is the wavelength, 

- R is the distance to the far field, 

- 0 is the half angle or cone angle from the axis to q. 

The hexagonal array is not quite circular, but in the limit it 
approaches a circle. In this case the constant 1.22 is simply 

replaced by F, a constant to be determined by experiment. 

By using Brigham's text [Ref. 19] it is possible to convert 
equation (7) into the form: 

(8) 0 = F ARCSIN (Ln/NAT) where, 
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- 0 is the angular measure of the side lobe from the main lobe, 

- L is the wavelength of interest = 6.86 cm, 

- NAT is the total width of the sampled grid = (ND/X), 

- N is the width of the input window (grid units;1 28,256 or 512), 

- D is the actual dimension of the physical array (cm), 

- X is the dimension of the input array in grid units, 

- n is the output number of grid units to the side lobe of 

interest=(rz/w), 

- r is the distance to the output side lobe of interest (cm), 

{NOTE: X and D must correspond; r and n must correspond}, 

- z is the width of the output window in grid units (< N/2), 

- w is the width of the input and output windows in centimeters. 

As an example, consider Figure 23. For this example, D = 
38.735 cm (measured between the outer edges of the farthest 2 
speakers), w = 10.3 cm, N = 256 grid points and z = 128 grid points. 
To calculate the angular extent of the main lobe in real space, first 
measure the distance (in cm) from the corner of the output quadrant 
to the edge of the main lobe, along the horizontal axis (Fig. 23). 
Here, r = 0.7 cm and X = 40 grid points. Therefore, equation (8) 
becomes: 

(6.86)(0.7)(1 28/10.3) cm 

(9) 0 = F x ARCSIN = F x 13.9° . 

(256)(38.735)/40 cm 
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Figure 23. Output Quadrant of Hexagonal FFT 

The half-angle of the actual array is only 12° as measured on 
Figure 8. This means that F * 0.86 or 1/F *1.16 which is very close 
to 1.22. Using the value F = 1.16, then the 90° point can be 
calculated for the output field. In this case, 

(256)(38.735) 

(10) Sin {(1 .1 6)(90°)} = r = 2.81 cm, 

(6.86)(1 28/1 0.3)(40) 
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which means everything outside the corner radius of 2.81 cm can be 
ignoredas it is non physical . 

The output FFT's were correlated to the actual lobe patterns of 
Figures 7 and 8, the apex and side mounted patterns for the 
hexagonal array, and the results are plotted in Figures 24 and 25 
and tabulated in Table 1. The amplitudes were normalized by 
assigning a value of unity to the main lobe; the angles are measured 
counter-clockwise from the standard X axis; the amplitude scale is 
in decibel units. 

The measured amplitude was calculated using INVLOG (dB/20), 
where 'dB' is the normalized decibel level estimated from Figures 7 
and 8. The computed amplitudes relied on the formula 20 LOG 
"Q/Max" [see Appendix B, lines 1390 and 2380 for definitions of Q 
and Max], 

The angles were computed using a simpler form of equation (8); 

(11) 0 = ArcSin { (K-1) L/NAt }, where K is an index from 

1 to N/2, and essentially allows the computation of the angle at 
each grid point chosen. The measured angles were taken directly 
from Figures 7 and 8. 
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Fig. 24. Polar Plot of Apex Mounted Hexagonal Array 
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Fig. 25. Polar Plot of Side Mounted Hexagonal Array 
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TABLE 1. COMPARISONS BETWEEN ACTUAL LOBE PATTERNS AND 
COMPUTED LOBE PATTERNS. 

APEX MOUNT ( Fig. 7 vs. 24 ) 

Measured Computed 





0 


amplitude 


0 


amolitude 


1 st min 


1 1 ° 


.025 


12° 


.002 


1 st max 


1 6° 


.126 


1 7° 


.128 


2 n d min 


21° 


.056 


23° 


.007 


2 n d max 


26° 


.095 


28° 


.108 


3 r< ^ min 


33° 


>.01 


38° 


.006 


3 r d max 


42 ° 


.056 


41 0 


.027 


min 


51 0 


.018 


46° 


.007 


4 th max 


64° 


.056 


62° 


.148 






SIDE MOUNT 


( Fig. 8 vs. 25 ) 




Measured 




Computed 




0 


amolitude 


0 


amplitude 


1 st min 


12° 


.050 


14° 


.048 


1 st max 


18° 


.158 


19° 


.151 


2 nc * min 


27 ° 


>.01 


28° 


.01 


2 n d max 


35° 


.039 


36° 


.030 


3 f d min 


52° 


>.01 


46° 


.024 


3 r< ^ max 


61 0 


.025 


56° 


.029 
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The side mounted section corresponded to measurements along 
the X axis, and the apex to measurements along the Y axis. The error 
of the measurements was grossly estimated as approximately 10% 
for the minimums, as they tended to be obscured, and 3% or less for 
the maximums. Some values exceed this amount considerably and 
require further investigation. 

Generally, however, the tabulated values agree strongly and 
suggest that this technique has considerable credibility at 
predicting side lobe structures of untested arrays. 
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VI. CONCLUSIONS 



This thesis demonstrated that increasing acoustic planar array 
speaker density from 25 speakers/1452 cm 2 (1.72 x 10' 2 
speakers/cm 2 ) to 19 speakers/943 cm 2 (2.02 x 10' 2 speakers/cm 2 ) 
effectively suppressed the lobe spreading from a 70° cone angle to 
less than a 45° cone angle, with a reduction in side lobe intensity of 
between 4 and 10 dB. 

Graphical profiles demonstrated the great variability and 
volatility of turbulence structures in the atmosphere and underline 
the need for better tools to measure and understand this structure. 
These profiles identify gravity waves, thermal plumes, boundary 
layers and other temperature induced structures. Particulary 
important, this device allowed tentative quantification of altitude, 
thickness and relative intensity of mountain boundary layers over a 
potential DOD laser site. 

Independent verification of the C T 2 values is important and is 

currently under study. It will further enhance the utility of this 
device. Due to the strong exponential dependence of acoustic 
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attenuation (3) with increasing frequency, future improvements 
should also include lower frequency drivers (2-3 kHz) of similarly 
high efficiency, to improve the range of the device. 

It was shown that acoustic array beam patterns could be 
successfully modeled on a computer by means of a Fast Fourier 
Transform algorithm and tested for optimal design. Future 
improvements to this technique include a method to stochastically 
simulate turbulence and model the beam structure as it passes 
through this artificial turbulence. Calculation of main lobe energy 
versus total energy for various arrays should also be studied. 
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APPENDIX A 



List of Material Sources 



Disclaimer: These sources are listed only in case exact 

duplication of this thesis is required at some future date. Other 
equally competitive sources for each item do exist and their 
exclusion is not a reflection of product quality in any way. 



Barium impregnated foam: 

1 . Barrier Corp. 

Tigard, Oregon 
(503) 639-4192 
POC: Mark Dove 

Cylindrical Enclosure: 

1 . Rubbermaid, Inc. 

1 147 Akron Rd 
Wooster, Ohio 44691 
(216) 264-6464 

Speakers: 

Motorola Inc. 

4800 Alameda Blvd N.E. 
Albuquerque, N. M. 87113 
(505) 822-8801 



2. Latchford Corp. 

(714) 734-3660 
POC: Tom Pelligreeni 



2. Traffic Products Co. 
Citrus Heights, Ca 
(916) 723-2958 
POC: Michael Fish 



Temperature/humidity Indicator 
Weathermeasure 
P.O. Box 41039 
Sacramento, CA 95841 
(916) 481-7565 
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APPENDIX B 



Fast Fourier Transform Program 



! Program plots either the input points, the FFT or the 
Polar Plots of a 19 element hexagon array or a 25 
! element square array. 

! This program is interactive. 

1 dlw Dec. 18, 86/LRM modified in OCT 87 
COM Re( 1:512) ,Im( 1:512) 

INTEGER M.N 

DIM Real(512 ,512 ) ,Imag(512 ,512 > 

REAL Amp( 1 00 ) ,Ang( 1 00 ) 

! 

PRINT "ENTER THE GRID SIZE N ; N= 512,256 OR 128" 

INPUT N 
PRINT 

i 

N2=N/2 

M = LOG( N )/L0G( 2 ) 



PRINT 


"ENTER 


’ 1 ’ 


TO 


SEE 


THE 19 


ELEMENT 


HEXAGON. “ 


PRINT 


"ENTER 


1 n » 

L. 


TO 


SEE 


THE 25 


ELEMENT 


SQUARE . " 


INPUT 

PRINT 

PRINT 


X 
















PRINT 


"ENTER 


’ 1 ’ 


TO 


SEE 


A 


PLOT 


OF THE 


INPUT POINTS 


PRINT 


"ENTER 


’ 2 1 


TO 


SEE 


A 


PLOT 


OF THE 


FFT" 


PRINT 


"ENTER 


» “7 1 
vJ 


TO 


SEE 


A 


POLAR 


PLOT OF INTENSITY" 



INPUT A 

PRINT 

PRINT 

PRINT "ENTER THE GRID SIZE OF EACH ELEMENT, R." 
PRINT " USUALLY R=4 OR 8." 

INPUT R 

PRINT 

PRINT 

PRINT "ENTER ’7’ TO SEE THE FIRST QUADRANT ONLY" 
PRINT “ ELSE ENTER ’0’ . " 

INPUT C 

i B= THE SIZE OF THE OUTPUT WINDOW 
B=N2/2 

I 

IF X= 1 THEN 

i 

F = 2 . 1 *R*C05( .5235988 ) 

FOR I=N2-6*R TO N2+6*R 
FOR J=N2-6*R TO N2+6*R 

IF ( ( I+4»R-N2 ) A 2 + ( J-N2 ) A 2 ><R*R THEN 
Real< I , J >= 1 .000 
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460 

470 

480 

480 

500 

510 

520 

530 

540 

550 

550 

570 

530 

550 

600 

610 

620 

630 

6 40 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 



FOR K= I TO 4 

Real( I+2*K*R , J )=Real< I ,J ) 

Real( I+2*K*R-R ,J+F >=Real< I , J ) 

Real< I + 2*K«R-R , J-F > = Real< I ,J ) 

NEXT K 

FOR K-l TO 3 

Real< I+2*K*R ,J-2*F )=Real< I , J ) 

Real! I +2 *K *R , J+2 *F )=Rea 1 < I , J ) 

NEXT K 
END IF 
NEXT J 
NEXT I 

I 

END IF 

i 

IF X=2 THEN 

i 

K - N 2 + 5 * R 
L=N2-5*R 
FOR 1=1 TO N 
FOR J=l TO N 

IF I<K AND I >L AND J <K AND J>L THEN 
Cy=N2-4*R 
FOR P-l TO 5 
Cx=N2-4«R 
FOR 0=1 TO 5 

IF I<Cx+R AND I;Cx-R AND J>Cy-R AND JxCy+R THEN 
IF R*R><( I-Cx>*2 + < J-Cy) A 2) THEN 
RealfJ.I >=1 .000 
END IF 
END IF 
Cx=Cx +2 *R 
NEXT 0 
Cy=Cy + 2 *R 
NEXT P 
END IF 
NEXT J 
NEXT I 

i 

END IF 

PRINTER IS 701 
! 

! This section plots the input points only 
! 

IF A= I THEN GOTO 2240 

i 

! THIS SECTION IS THE FAST FOljRIER TRANSFORM 

l 

T I =T I MEDATE 
FOR 1=1 TO N 
FOR J=1 TO N 
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970 

980 

990 

I 000 

1010 

I 020 

1 030 

1 040 

1 050 

1060 

1 070 

1 080 

I 090 

1 100 

1110 

1 120 

1 1 30 

1 1 40 

I 150 

I 1 60 

1 1 70 

I 1 80 

I 1 90 

1 200 

1210 

1220 

I 230 

I 240 

1 250 

I 260 

1 270 

1230 

1 290 

1 300 

1310 

1 320 

1 000 

1 340 

1 350 

1 360 

1 370 

1 380 

1 390 

I 400 

1410 

I 4 20 

1 430 

1 440 

I 450 

1 460 

I 470 



Re ( J ) = Rea 1 ( I ,J ) 

Im( J )=Imag( I , J ) 

NEXT J 

i 

CALL Fft(M) 

i 

FOR J=1 TO N 

Real( I , J )=Re( J ) 

I m a g ( I ,J ) = I m ( J ) 

NEXT J 
NEXT I 

i 

FOR J=1 TO N 
FOR 1=1 TO N 

Ret I )=Rea 1 ( I ,J ) 

In ( I )=Imag( I , J ) 

NEXT I 

i 

CALL Fft(M) 

i 

FOR 1=1 TO N 

Real( I , J >=Re( I ) 

I mag( I , J ) = Im( I ) 

NEXT I 
NEXT J 
T2=TIMEDATE 

! PRINT I NT ( 1 00( T2-T1 >/60 )/ 1 00 , "MINUTES" 

l 

IF A=2 THEN GOTO I960 
! 

! This section graphs the polar plot of the FFT. 

I 

61 =N2/2 
62=1 

1 0= The actual dimension of the array(cm) 

D=38 . 735 

Pma>. = ( 1 +INT< D*N/(68. 6*R ) ) ) 

i 

! The ’ K 1 loop does the side and apex mounts 
! 

R 1 =Rea 1 ( 1,1) 

Ia= Imag( 1,1) 

lia.\ = SQR(Rl # Rl + Ia*Ia ) 

i 

FOR K= 1 TO 2 

i 

IF X=1 AND K= 1 THEN 

PRINT " HEXAGONAL. ARRAY , SIDE MOUNTED 

END IF 

IF X= 1 AND K=2 THEN 

PRINT " HEXAGONAL ARRAY, APEX MOUNTED 
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I 480 

1490 

I 500 

1510 

1520 

1 530 

1 540 

I 550 

I 5G0 

1 570 

1580 

1590 

I 600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1 730 

1740 

1 750 

1 760 

1 770 

1 780 

I 790 

1 800 

1810 

I 820 

I 830 

I 840 

1850 

1 860 

1 870 

1880 

1890 

I 900 

1910 

1 920 

1930 

1 940 

1 950 

1 960 

1 970 

I 980 



END IF 

IF X = 2 AND K= 1 THEN 

PRINT “ SQUARE ARRAY, SIDE MOUNTED 

END IF 
! 

1 This section computes the ape* square array 
i 



IF X=2 AND K=2 THEN 

PRINT * SQUARE ARRAY, APEX MOUNTED" 

1 The 'L’loop does only the apex mounted square 
FOR L= 1 TO N2 
R 1 =Rea 1 ( L ,L ) 

G=Imag( L ,L ) 

Q=SQR< R1 *R1+G*G ) 

IF L > P m a a THEN GOTO 1660 
Amp ( L )=Q/Max 

Ang( L )=ASN< < (L-l ) *68 . 6 * R >/< D*N ) ) ♦ 1 80/P I 
! PRINT L ,Q/Max ,ASN( < (L-1 >*68. 6*R >/<D*N ) >*180/3. 1415927 
NEXT L 

CALL Polarplt(Amp< * ) ,Ang( * ) ,Pmax ) 

IF X =2 THEN GOTO 2590 
END IF 



i 



FOR 1=1 TO B 1 
FOR J=1 TO B2 
IF B 1 = 1 THEN P=J 
IF B2= 1 THEN P=I 
R 1 =Rea 1 ( I ,J ) 

6= I mag( I , J ) 

Q = SQR< R1 *R 1 +G*G ) 

IF P>Pmax THEN GOTO 1840 
Amp( P )=Q/Max 

Ang( P )=ASN( < (P-1 ) *68 . 6 *R >/ ( D *N > ) * 1 80/P I 
! PRINT P, Q/Max ,ASN(((P-1 > *68 . 6 *R > / ( D*N > > « 1 80/3 . 1415927 
NEXT J 
NEXT I 

I 

CALL Polarplt(Amp( * > ,Ang( * ) ,Pmax ) 

! 

B2=B1 
B 1 = 1 
NEXT K 
! 

IF A = 3 'THEN GOTO 2590 
! 

! This section centers the FFT onto the output grid. 

i 

IF C=7 THEN ! USE ONLY TO GET ONE QUADRANT OF THE GRID 
GCLEAR 
GINIT 
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1 330 
2000 
2010 
2020 
2030 
2040 
2050 
2050 
2070 
2080 
2090 
2100 
2110 

2 I 20 
2 I 30 
2140 
2150 
2150 
2170 
2180 
2 1 90 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2 j o 0 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 



GRAPHICS ON 
0 1 EUIPORT 30,130,0,100 
WINDOW 1 ,N2 , 1 ,N2 
61 =N2-1 
62=0 

GOTO 2340 
END IF 
! 

1 This section centers the output, 
t 

FOR 1=1 TO N2-I 
FOR J=1 TO N2-1 
R I = Rea 1 ( I , J ) 

I a = I mag( I , J ) 

Real ( N2+I-1 .N2 + J-1 > = R1 
Real ( N2- 1 ,N2 + J-1 ) = R1 
Real ( N2-I ,N2-J )=R1 
Real ( N2 + 1 - I ,N2-J ) = R1 
Imag(N2-I ,N2-J )=Ia 
Imag( N2 + I- 1 ,N2-J )=Ia 
Inag( N2+I-1 .N2+J-1 )=Ia 
I rnag( N2- 1 , N2 + J- 1 )=Ia 
NEXT J 
NEXT I 

l 

61=6 

62=6 

! 

GCLEAR 

GINIT 

GRAPHICS ON 
! 

VIEWPORT 30,130,0,100 
WINDOW N2-B ,N2 + 6 ,N2-B ,N2 + B 

l 

FOR I =N2-B 1 TO N2+B2 
FOR J =N2-B 1 TO N2 + B2 
R 1 =Rea 1 ( I ,J ) 

Ia=Inag( I , J ) 

Q = SQR(R1*R1 + Ia*la ) 

IF A= 1 THEN 

Z=( 12*0+1 >/1G 1 USE ONLY FOR PLOTTING INPUT POINTS 

END IF 
IF A=2 THEN 

Z = 1+LGT(Q/Ma.O/D ! ONLY FOR DECADE INTENSITY SPREAD 
END IF 

IF 2 <0 THEN Z = 0 
IF Z>1 THEN Z=l 
AREA INTENSITY Z ,Z ,Z 
MOVE I ,J 

RECTANGLE 1 , 1 .FILL 
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2500 

2510 

2520 

2530 

2540 

2550 

2550 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 



NEXT J 
NEXT I 

i 

PRINT "The size of the output window is",B,"X",B 
PRINT "The size of the input window is " ,N , “X " ,N 
PRINT "The radius of each element is " ,R 

j 

1 Pause statements allow insertion of parameters w/o 
! redoing the FFT calculations. 

PRINTER IS CRT 

PAUSE 

PAUSE 

END 

i 

1 This section is the FFT subroutine (from INFOTEK). 
i 

SUB Fft( INTEGER M> 

CON Re( 1:512) ,Im( 1:512) 

! FFT SUBROUTINE FRON INFOTEK 29 AUG, 1985: 

! USE COMMON FOR FASTER EXECUTION 3.09 sec vs 3.59 

INTEGER I , J ,K ,L ,N , Ip ,Le ,Lel 

N=2~M 

J=1 

FOR 1=1 TO N-l 
IF I<J THEN 
T=Re( J ) 

Re ( J >=Re< I ) 

Re( I ) = T 
T=Im( J ) 

Im( J ) = Im( I ) 

Im( I ) = T 
END IF 
K=N DIU 2 
WHILE K< J 
J = J-K 
K-K DIU 2 
END WHILE 
J = J+K 
NEXT I 
Le= 1 

FOR L= 1 TO M 
Le 1 =Le 
Le=Le+Le 
Ure= 1 
U l m=0 
Ang=P I /Le 1 
Wre=C0S ( Ang ) 

W l m = S I N( Ang ) 

FOR J=1 TO Le! 

FOR I=J TO N STEP Le 
Ip=I+Le I 
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3010 

5020 

3000 

3040 

3050 

3080 

3070 

3030 

3090 

3100 

3110 

3120 

3130 

3U3 

3150 

3160 

3170 

3 I 80 

3190 

3200 

32 10 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

331 0 

3320 

3330 

0 0 4 0 

oo 50 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

34 50 

3460 

3470 

3480 

3490 

3500 

3510 



Tre = Re( Ip ) »Ure-I m( Ip ) «Uin 
T im = Re( Ip ) ♦Uin+ Im( Ip )*Ure 
Re( Ip ) = Re< I )-Tre 
In< Ip ) = Im( I )-T in 
Re< I )=Re( I ) + Tre 
I m( I >= In< I ) + T in 
NEXT I 

T=Ure*Wre-Uin*Wim 
Uim=Ure*Win+Uin«Wre 
Ure = T 
NEXT J 
NEXT L 

Npscale*! . /N 
MAT Re= <Npscale)*Re 
MAT In= ( Np scale )*In 
SUEEND 
! 

1 This section is the polar plot subroutine written 
! by Paul Davi son . 



SUB Polar p 1 1 ( Amp ( * ) ,Ang( * ) , Pmax ) 

REM 

REM THIS MODULE MAPS THE CALCULATED, NORMALIZED 

REM DSCS ARRAY DS< * ) ON LOG NORMAL POLAR PLOTS. 

REM 

i 

i ++++ INITIALIZING THE PLOTTER 

i 



GINIT 

PLOTTER IS CRT , " INTERNAL 

GRAPHICS ON 

GCLEAR 

i ++++ SCALING THE GRAPH 

VIEWPORT 25 , 100 ,20 ,95 
SHOD -40 ,40 ,-40 ,40 



LINE TYPE 4 
R l ng= 1 0 

AXES Ring .Ring ,0 ,0 ,5 ,5 ,2 

LINE TYPE 3 

DEG 

LORG 7 
CSIZE 3 
FOR R=0 TO 4 
LINE TYPE 1 
MOVE 0 ,R*Ri ng 
LABEL R*Ri ng-40 
LINE TYPE 3 
FOR Ang 1 e=0 TO 360 

PLOT R»Ring«COS( Angle ) ,R»Ring*S INI Angle ) 
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352a 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

j6 j0 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 



NEXT Angle 
PENUP 
NEXT R 

FOR Ang 1 e=30 TO 150 STEP 30 
IF Angle<>90 THEN 

PLOT 40*COS( Angle ) , 40 *S IN( Angle ) ,-2 
DRAW 40*COS( 1 80+Angl e ) ,40*SIN( 1 80+Angle ) 
END IF 
NEXT Angle 
PENUP 

i 

LIME TYFE I 
DEG 

FOR F=Pmax TO 1 STEP -1 
R=40+20*LGT(Anp(P ) ) 

IF R>=0 THEN 

PLOT R*COS(Ang(P ) ) ,R*S IN( Ang< P ) ) 

END IF 
NEXT P 

FOR P=l TO Pmax 

R=40+20*LGT(Amp(P ) ) 

IF R> = 0 THEN 

PLOT R*COS( Ang( P ) ) ,-R*SIN< Ang< P ) ) 

END IF 
NEXT P 
! 

CLIP OFF 
! 

LORG 2 
MODE 40,0 
LABEL " X" 

LORG 8 
MODE -40,0 
LABEL "-X " 

LORG 6 
MODE 0,-40 
LABEL “-Y" 

i 

CLIP ON 
! 

DUMP GRAPHICS #701 
GRAPHICS OFF 
PRINT CHR$( 12 ) 

SUBEND 
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